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Abstract 

We calculate the mass and decay constant of / = ^ scalar mesons 
composed of quark-antiquark pairs based on QCD sum rule. The quauk- 
antiquark pairs can be sq or qs {q = u, d) in quark model, the quantum 
numbers of spin and orbital angular momentum are S = 1, L = 1. We 
obtain the mass of the ground sate in this channel is 1.410 it 0.049GeV. 
This result favors that Kg (1430) is the lowest scalar state of sq or qs. 
We also predict the first excited scalar resonance of sq is larger than 2.0 
GeV. 

1. Introduction 

Glueball and scalar mesons should exist according to QCD and quark model. 

Some scalar mesons below 2 GeV have been observed, such as, i) for isospin 

/ = 0,1 states: /o(600) or a, ao(980), /o(980), /o(1370), /o(1500), /o(1710); 

ii) for / = I states: fi:(900) and K*{U30) [1 H El i]. The number of these 

scalar mesons exceeds the particle states which can be accommodated in one 

nonet in the quark model. It is believed that there are two nonets below and 

above 1 GeV |H1 E] • The components of the meson states in each nonet have 

not been completely determined yet. For the scalar mesons below 1 GeV there 

are several interpretations. They are interpreted as meson-meson molecular 

states |7j or multi-quark states qqqq |H1, etc.. However, from the theoretical 
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point of view there must be quark-antiquark SU(3) scalar nonet. Therefore it 
is important to determine the masses of the ground states of qq with quantum 
number = 0^ based on QCD. For isospin J = 0, 1 states different quark 
flavor may mix, and scalar qq states may also mix with scalar glueball if they 
have the same quantum number of J^*" and similar masses 0-^3- Some 
authors have tried to determine the mixing angles of the glueball with qq 
scalar mesons by using decay patterns of some scalar mesons These 
works imply that glueball possibly mix with qq scalar mesons. For J = 1/2 
states, they cannot mix with glueball because they have strange quantum 
number. The physical state is directly the sq and qs bound state. Therefore 
the mass of the ground state of sq or qs can be determined without necessity 
for considering mixing effect. 

QCD sum rule is a powerful tool to calculate hadronic nonpertubative 
parameters based on QCD j^. It has been used to calculate the masses and 
decay constants of 0^^,1^^,2++ mesons before and give satisfactory results 

dHllinillllEI!. 

In this paper, we calculate the mass and decay constant of / = 1/2 scalar 
meson with QCD sum rule. We find that it is impossible to obtain sq scalar 
meson mass below 1 GeV from QCD sum rule. The most favorable result for 
the mass of sq scalar meson is 1.410±0.049GeV. Therefore, if k(900) is sq scalar 
bound state, this would be a big problem for QCD. This problem can be solved 
by assuming that k(900) is irrelevant to sq scalar channel, < 0|sg|/t(900) >~ 0, 
and -^"0(1430) is the scalar ground state of sq or qs. With this assumption, 
calculation based on QCD will be consistent with experiment. Therefore, our 
result favors that i^'Q(1430) is the lowest scalar bound state of sq. 

If this is correct, then from the approximate SU(3) flavor symmetry, the 
masses of the other = 0^ mesons in the scalar nonet should be slightly above 
or below 1.4 GeV. This result would imply that scalars with masses below 1 
GeV are not dominated by quark-antiquark pairs. This is consistent with 
the calculation of lattice QCD which implies that a nonet of quark-qntiquark 
scalars is in the range 1.2-1.6 GeV |22j . 

The remaining part of this paper is organized as follows. In Section 2, we 
briefly introduce the process to calculate the scalar meson with QCD sum rule 
and get the Wilson coefficients for the corresponding two-point scalar current 
correlation function. Section 3 is devoted to numerical analysis and conclusion. 
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2. The method 

To calculate the mass of scalar sq or qs meson, the two-point correlation func- 
tion should be taken as 

n(g2) = , j rfW^-(o|r{j(x)j+(o)}|o) (1) 

where = s(x)g(x), j^(0) = g(0)s(0). 

On one hand, the correlation function can be expressed based on the dis- 
persion relation in terms of hadron states 

irdjhji^ (2) 
-K J s — q^ 

where Im^{s) is the imaginary part of the two-point correlation function, which 
can be obtained by inserting a complete set of quantum states J2 |'^)('^| into 
Eq.((H). The result is 

2/^n(s) = 5:27r5(s-m2)(0|j(0)|n)(n|j+(0)|0) (3) 

n 

For the scalar states S, its decay constant fs can be defined through 

{0\j{0)\S) = msfs (4) 

where m5 is the mass of the scalar state. Based on Eq.©- Eq.(@]), and ex- 
plicitly separating out the lowest scalar state, the correlation function can be 
expressed as 

li\q^) = + - / 5 

rrig — q'^ ti Js° s — q^ 

where p^{s) expresses the contribution of higher resonances and continuum 

state, s° is the threshold of higher resonances and continuum state. 

On the other hand, the correlation function can be expanded in terms of 

operator-product expansion at large negative value of q^. 

n«^^(g2) = z|d^a;e^'^-(0|T{j(a:)j+(0)}|0) 

+ C6(o|^'r$fr'$|o) + --- (6) 

where Cj, i = 0, 3, 4, 5, 6, ■ ■ ■ are Wilson coefficients, I is the unit operator, 
^f^' is the local Fermion field operator of light quarks, C^p is gluon strength 
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tensor, F and F' are the matrices appearing in the procedure of calculating the 
Wilson coefficients. 

For convenience later, we reexpress the above equation as 



s — 



(7) 



where Ps""^, ■ ■ ■ , Pe""^? ' ' " ^-^^ contributions of condensates of dimension 3,4, 
5, 6,--- in Eq.®. 

Matching Il^{q^) with W^^{q^) we can get the equation which relates mass 
of scalar meson with QCD parameters and a few condensate parameters. In 
order to suppress the contribution of higher resonances and that of condensate 
terms, we make Borel transformation over in both sides of the equation, the 
Borel transformation is defined as 

2. i-Q'^T 5" 



5, 



lim 



n — > oo 
—>■ — oo 

2 _ nj2 



[n — 1)! d{q 



2\ri 



-q^/n = 

After assuming quark-hadron duality, i.e., by assuming that the contribution 
of higher resonance and continuum states can be approximately cancelled by 
the perturbative integration over the threshold s° [21], the resulted sum rules 
for the mass and decay constant of the scalar meson are 



ms 
fs 




1 

ms 



where 
Ri = 



(8) 
(9) 



R2 



,0 



(9M2 



(10) 



s)e 



71 7(mi+m2)2 

+M25pr"^ + M'^Bpl'"'^ 



+ M'Bp'l"''^ + M'Bpl 



nonp 



where Bp^°^^ express Borel transformation of pi°^^ ■, M is Borel parameter, and 
mi and 1712 are the masses of the two light quarks. 
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We need to calculate the Wilson coefficients in Eq.® to get the mass and 
decay constant of scalar meson. Collecting the contribution of diagrams in 
Fig.l, we get the result of Bp^°"'^ which is listed in Appendix. 




Figure 1: Diagrams for the contribution to Wilson coefficients, (a): diagrams con- 
tribute to unit operator; (b): diagrams contribute to bi-quark operators 'l&(a;)^'(0); 
(c): diagrams contribute to G^^G""^"; (d): diagrams contribute to quark- gluon mix- 
ing 'I'(x)^'(0)G^;^; (e): diagrams contribute to four-quark operators (^^f)^ 

3. Numerical analysis and conclusion 

The numerical parameters used in this paper are taken as jTHl |^ 
{qq) = -(0.24 ± 0.01GeV)3, {ss) = ml{qq) 
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as{GG) = 0.038GeV^ ^(^aT^) = m2(^^) 



= 6.0 X 10-^GeV^ = 0.8 ± 0.2GeV^ 
m, = 0.14GeV, m„ ^ = 0.005GeV (12) 

For the choice of Borel parameter M^, as in [HI 1211, define fthcorriM"^) as 
m(M^) in Eq.® without the continuum contribution (s° = oo) and mnopower{M'^) 
as m{M'^) in Eq.® without power corrections, then define fnopower{M'^) as 
m(M2)/m„opo^er(M2) and fcont as m{M'^) / fthcorr{M'^). To get rehable predic- 
tion of the mass in QCD sum rule, fcont should be limited to above 90% to sup- 
press the contribution of higher resonance and continuum, and fnopoweriM"^) 
be limited to less than 10% deviation from 1, which can ensure condensate 
contribution much less than perturbative contribution. 

There are two low mass scalar meson states with isospin 1 = 1/2 and 
strange number 15*1 = 1 found in experiment. They are fi:(900) with mass 
m« about 800 ~ QOOMeV 013111, and K*{U30) with mass m{K*{U30)) = 
1.412 ± 0.006GeV [I]. In theory, taking appropriate value for the threshold 
parameter s", one can separate out the contribution of the lowest resonance in 
QCD sum rule. We vary the value of the threshold parameter s°, and find that 
it is impossible to obtain the mass of k(900) with the sum rule in Eq.®. There 
is no stable 'window' for the Borel parameter in this mass region. Therefore, if 
fi:(900) is the lowest scalar state in the sq channel, it would be a big problem for 
QCD sum rule. However, if we increase the value of i.e., for s° = 4.0 ~ 4.8 
GeV^, we does find the stable 'window' for Borel parameter, which is shown 
in Fig 121 The resulted stable window is in the range 1.0 < < 1.2 GeV^. 

Fig|2fa) shows that between the arrows A and B, both the contributions 
of condensate and higher resonance are less than 10%. So in this region, 
the operator product expansion is effective, and the assumption of quark- 
hadron duality does not seriously affect the numerical result, which means 
that QCD sum rule can give reliable prediction in this parameter space. For 
s° = 4.0 ~ 4.8 GeV^, the mass of scalar sq meson in QCD sum rule is 

m{sq) = 1.410 ± 0.049 GeV (13) 

where the error bar is estimated by the variation of Borel parameter in the 
range 1.0 < < 1.2 GeV^, the variation of s° within 4.0 ~ 4.8 GeV^ 
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Figure 2: (a) The region between the arrow A and B is reliable for determining the 
mass { = AAGeV"^). (b) The curves correspond to the mass of scalar meson for 
the continuum threshold s° = 4.0GeF^,s° = AAGeV"^ , A.SGeV'^ , respectively. The 
central one is for = AAGeV'^. 



the uncertainty of higher «<, correction for the perturbative diagram and the 
condensate parameters. The variation of Borel parameter yields ±1.8% uncer- 
tainty for the mass, s° yields ±2.0%, as correction gives ±2.2%, the uncertainty 
caused by the condensate parameters is less than 0.6%. All the uncertainties 
are added quadratically. 

The energy scale for the as{fi) correction is taken to be = M. In the 
stable window, the range of Borel parameter is 1.0 < < 1.2GeV^, therefore 
as{M) ~ 0.5. We checked that the contribution of the correction at first 
order is about 2.2%, which is not large. This can be understood because 
most contribution of the correction is cancelled between the numerator and 
denominator of Eq.®. We use 2.2% to estimate the uncertainty caused by 
the higher order corrections. 

On one hand, it is impossible to obtain the mass of lower scalar state /€(900) 
from QCD sum rule for sq channel. If regard k(900) as sq scalar bound state, 
it would be a big problem for QCD. On the other hand, QCD sum rule can give 
most favorable mass which is consistent with the mass of Kq{1430). Therefore 
it is acceptable to assume that fi:(900) is irrelevant to sq scalar bound state, 
and 

< 0|sg|K(900) >~ (14) 
With this assumption, Kq{1A30) can be accepted as the lowest scalar bound 
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state of sq. Then there will be no problem between QCD and experiment. 

One may still be afraid that there are contributions of the lower mass state 
k(900) mixed in the result of eq. ()13|l in fact. If this is indeed the case, the 
result of the sum rule may be some weighted average of the two resonances 
of fi;(900) and i^'o(1430). Therefore this situation should be carefully checked. 
Because the sum rule for the mass of the scalar bound state in eqs.dHI), (fTHj) and 
(fTT|) includes the spectrum integration j^^^j^^^y ds, in principle one can lower 
the value of to separate the lowest bound state. Therefore, we checked what 
result for the mass can be got by lower the value of s° within the stable window 
1.0 < < 1.2GeV2 selected in FigEfa). The resuh is shown in FigEl It 
shows that for any value of the possible mass is large than 960MeV, 

m{sq) > 960 MeV (15) 

Therefore the possible effect of k(900) can be safely ruled out in the sum rule 
result in eq. p3p . Note that the most recent experimental result for the mass 
of fi;(900) from E791 collaboration is m« = 797 ± 19 ± 42MeV |SI. 

m(GeV) 
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1 . 5 

1 . 4 

1 . 3 

1 . 2 

1 . 1 



12 3 4 
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Figure 3: The possible mass result by varying the value of the threshold s^. The 
solid curve is for the Borel parameter = l.OGeV^, and the dashed one for 
M2 = 1.2GeV2. 

If i^'Q(1430) is the ground state of sq or qs, from the approximate SU(3) 
flavor symmetry, the masses of the other = mesons in the scalar nonet 
should be also around 1.4 GeV. This implies that the scalars with masses 
less than 1 GeV, i.e., /o(600), ao(980), /o(980) etc., can not be dominated by 
quark- antiquark bound states. This is consistent with the calculation of lattice 
QCD which implies that a nonet of quark- antiquark scalars is in the region 
1.2-1.6 GeV [221- 
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Our result can be further checked by experiment. From the threshold 
parameter s", we can predict that the mass of the first excited resonance in sq 
scalar channel should be larger than \f^, that is 

m*(K*) > 2.0 GeV (16) 

This prediction can be tested by experiment. 

Next we discuss the decay constant of the two-quark scalar bound state 
sq. From the above analysis, we take the threshold parameter s° = 4.0 ~ 4.8 
GeV^. Consider ir*(1430) as the only resonance below 2 GeV in the sq scalar 
channel, we can obtain the decay constant of frQ(1430) as a function of Borel 
parameter (see eq.dH}). The numerical result is shown in Fig|3 

/a-* (1430) (GeV) 

0.8 - 
0.6 - 

0.2 - 
0.6 0.8 1 1.2 1.4 1.6 1.8 

M^{GeV^) 

Figure 4: The decay constant of i('Q(1430) as a function of the Borel parameter 
Af2. The solid curve is for = 4.0GeV^, and the dashed one for = 4.8GeV^. 

Fig.4 shows that the decay constant is very stable. The determined stable 
'window' is still in 1.0 < < 1.2GeV^, where the continuum and condensate 
contribution are restricted to be less than 15% and 4%, respectively. Within 
this stable window, the decay constant of Kq{1A?)Q) is 

/(fs:*(1430)) = 427 ± 85 MeV (17) 

The variation of s° yields ±30% uncertainty for the decay constant, correc- 
tion gives ±20%, the uncertainties caused by the condensate parameters and 
the variation of Borel parameter are less than 0.3% and 0.1%, respectively. All 
the uncertainties are added quadratically to give the error bar in the above 
result. 

Again we should check what will happen if we consider two resonances 
/t(900) and 7^^0(1430) existing below 2 GeV in our sum rule analysis. Therefore 
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we add one more resonance into eq.©, then matching n'*(g^) with n'^*^^(g^) in 
eq. ((Z)) . By assuming quark-hadron duahty to cancel the contribution of higher 
resonance and continuum above 2 GeV, and making Borel transformation in 
both sides, we get the Borel improved matching equation 



m 



;i8) 



where R2 has been given in eq. (fTT|) . and msi, ms2 are fixed to be the masses 
of k{900) and i\:o*(1430), mgi = 900MeV, ms2 = 1410MeV. fsi and fs2 are 
the decay constants of the relevant scalar mesons. 

Differentiate both sides of eq.(fTH| with the operator d/dM'^, we can get 
another equation 



m 



where Ri is defined in eq-lfTUj). With eqs.(fTH|l and fT^. we can obtain 



fsi 



/; 



mli{ml2-ml^) 



e"'s2/ 



52 



:mi,2R2-Ri) 



ml^R2 - Ri] 



(19) 

(20) 
(21) 



From the above result we can perform the numerical analysis for the decay 
constants in the two- resonance ansatz. The numerical result is shown in Fig|31 
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Figure 5: The decay constants in two-resonance ansatz below 2GeV . The solid 
curve is for = 4.0GeV^, and the dashed one for = 4.8GeV^. (a) The decay 
constant of the low resonance k(900). (b) The decay constant of the higher resonance 

K^{im) 

From FiglHl we can see that both the two decay constants are unstable 
as a function of Borel parameter in the two-resonance ansatz. Adding the 
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lower resonance k(900) in the sum rule analysis for the sq channel spoils the 
stability existing in the one-resonance ansatz, which is shown in FigjH From 
the requirement of numerical stability of QCD sum rule, the numerical analysis 
of the decay constant does not favor to include k(900) in sq scalar channel. In 
addition, we can see from Fig|3^ that the decay constant of the lower scalar 
resonance fi:(900) tend to be zero at ~ 1.01 and 1.05 GeV. This is consistent 
with the requirement that < 0|sg|K(900) >~ in the one-resonance ansatz, 
where the stability window is located in the range 1.0 < < 1.2 GeV. 

Therefore, both the analyses of the mass and decay constant of sq scalar 
meson from QCD sum rule imply that fi;(900) is not dominated by quark- 
antiquark bound state, and the lowest sq scalar bound state is i^'Q(1430). The 
mass obtained from QCD sum rule is 

m(ifo(1430)) = 1.410 ± 0.049GeV (22) 

and the decay constant is 

/(ir*(1430)) = 427 ± 85MeV. (23) 

In summary, we calculate the mass and decay constant of scalar meson sq 
in QCD sum rule. Our result favors that -^"0(1430) is the ground state of sq 
scalar bound state. If this is correct, it would imply that scalar mesons below 
1 GeV are not dominated by quark-antiquark pairs. We also predict that the 
mass of the first excited resonance of sq scalar bound state is larger than 2.0 
GeV. 
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Appendix 

Borel transformed coefficients of perturbative and nonperturbative contribu- 
tions 13 pI"^^ in Eqs.firUj) and (fTTj) are listed below 

-3[(mi + - s\M-{mi - m-iY + s)(-(mi + m^Y + s) 
^ = ^ 8^^^ 
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Svr 3 vr ^ ^ ' 

where the term with as(yu) is the radiative correction to the perturbative con- 
tribution [20], and the scale is taken to be /i = M. 



Bpl°^P = ['iM^miml + ^M'^m\ml + mlml + ^M^{mi + 2m2)] 



e 



ml/M^ 



+?,M\2m, + m^Wd)^^ (25) 



Bpr' = 4vra,(GG){ -3K + ^2 



256e(Ki+"^2)'/A^2)M2mim27r2 

+ {3M'^mlm2 + SM^mlml + m\m\ + 3M'^(2mi + m^)) 
1 



288eKA^')M8m27r2 

+ {3M'^miml + 3M'^mlml + m\m\ + 3M^(mi + 2m2)) 
1 



288e(™i/^^')M8mi7r2 
—Vlm\{m\ — m2)^m2 + M^(— 7m^ + 26mim2 — 7m2) 



768e(('"i-'"2)VA^2)j^^4^^^27r 



2 



3(mi + ma)^ 



(mi+m2)2 ^ 128e[('^i+"^2)VM2]M47r2(mf + 2mim2 
m\m2{vn\ — m\m2 + vn^ — t)t^ 
+ 8g(t/M2)j^^47^2(^2 _ 2mxm2 + m| - t)2(mf + 2mim2 + m| - t) 

— {(mi — m2)^[4mi(mi — m2)^vrL2{vn\ — 2m\m2 + ^2 — t) 
+M^(3m^ — 16m'^m2 + 2^m\m^ — IGrriiml + ?)m\ — 
2>m\t + 14mim2t — 3m2t)]} 

\ } 

128e(('"i-'"2)'/*^')M67r2(mf - 2mim2 + - tf^ 

, ^ } (26) 

y^[(-(mi - m2)2 + t)(-(mi + m2)2 + t)\ 

-nonp ,= mi[-6M^ + m?m2 + 3M2mi(mi+m2)] 
Sps = ^(^^tT^')! i2eK/^^^)M8 

m2[— GM*^ + mimg + 3M2m2(mi + ma)] , 
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+36M^mi(mi + 2m2) + 84:M^ml{ml - mj) 
+15M'^mt{ml - ml)] ^ 



81e(™2/A^^)M8mi(-mf + mi) 
+[36M^m2(2mi + ms) + ml{ml - ml) 

-8AM\mt - mjml) - 15M\ml - mfm^)] 

81eK/^')M8m?(mf-mi)^ 
where mi = m^, m2 = m^. 
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